Quantum phase estimation with lossy interferometers 
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We give a detailed discussion ol optimal quantum states for optical two-mode interferometry in 
the presence of photon losses. We derive analytical formulae for the precision of phase estimation 
obtainable using quantum states of light with a definite photon number and prove that maximization 
of the precision is a convex optimization problem. The corresponding optimal precision, i.e. the 
lowest possible uncertainty, is shown to beat the standard quantum limit thus outperforming classical 
interferometry. Furthermore, we discuss more general inputs: states with indefinite photon number 
and states with photons distributed between distinguishable time bins. We prove that neither of 
these is helpful in improving phase estimation precision. 
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I. INTRODUCTION 

The strong sensitivity of certain quantum states to 
small variations of external parameters opens up great 
opportunities for devising high-precision measurements, 
e.g. of length and time, with unprecedented accuracy. A 
particularly important physical measurement technique 
is interferometry. Its numerous variations include Ram- 
sey spectroscopy in atomic physics, optical interferome- 
try in gravitational wave detectors, laser gyroscopes and 
optical imaging to name but a few. Understanding limits 
on its performance in realistic situations under given re- 
sources is therefore of fundamental importance to metrol- 
ogy. In this paper we examine the fundamental limits of 
the precision of optical interferometry in the presence of 
photon losses for quantum states of light with definite 
photon number. 

Optical interferometry aims to estimate the relative 
phase of two modes, or two "arms" , of the interferome- 
ter. This estimation process requires a certain amount of 
resources which is typically identified to be the number of 
photons, N , used for the measurement. The best preci- 
sion which can be obtained using classical states of light 
scales like the so-called standard quantum limit 

(SQL). Using non-classical states of light this precision 
can be greatly improved, ideally leading to Heisenberg- 
limited scaling, 1/7V [HQ- Indeed, recent years have seen 
many experimental proof-of-principle demonstrations of 
beating the SQL using quantum strategies in various in- 
terferometric setups [3, 0, H, H, 0, Unfortunately, 
highly non-classical states of light which potentially lead 
to Heisenberg-limited sensitivity are very fragile with 
respect to unwanted but unavoidable noise in experi- 
ments. In quantum-enhanced optical interferometry, the 
loss of photons is the most common and potentially the 
most devastating type of noise that one encounters. In 
particular it was noted that highly entangled quantum 
states, optimal for interferometry in the lossless case ~ 



NOON states [13] - are extremely fragile. Even for mod- 
erate losses they are outperformed by purely classical 
states [nl, [H, [ij, [13, d, In, 113 • A different approach 
has been taken in [9|, where the noise arising from im- 
perfect preparation of a state has been investigated. 

In [m , the first systematic approach was taken in order 
to determine the structure of optical states optimal for in- 
terferometry in the presence of losses. The best possible 
precision using input states with definite photon num- 
ber was given. In this paper we elaborate and extend the 
ideas presented in 15]. Our treatment is based on general 
quantum measurement theory [H, [l^ [13, [llj . The quan- 
tity of interest is the lowest possible uncertainty attain- 
able in parameter estimation, inversely proportional to 
the square root of the quantum Fisher information [13] ■ 
The quantum Fisher information depends only on the 
state of the system and not on the measurement proce- 
dure. We show that the optimization of the quantum 
Fisher information can be done effectively over the class 
of input states which have a definite photon number. 
Since these input states are subject to unavoidable pho- 
ton losses they will degrade into mixed states and their 
suitability for phase estimation is compromised. Our op- 
timization takes this into account yielding the most suit- 
able input states in the presence of photon losses lead- 
ing to the highest possible quantum Fisher information, 
and hence to the best possible precision. We give a de- 
tailed description of the noise model and calculate an 
analytic expression for the quantum Fisher information. 
The latter is shown to be a concave function on a convex 
set and therefore suited for efficient convex optimization 
methods. We numerically determine the optimal input 
states, compare them to alternative quantum and classi- 
cal strategies, and show that they can beat the SQL. We 
note that the corresponding precision, which lies between 
the SQL and the Heisenberg limit (depending on the loss 
rates) , defines the best possible precision for optical two- 
mode interferometry. In addition to this we discuss a 
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FIG. 1: Interferometric phase estimation scheme. An input 
state is fed into an interferometer consisting of two channels. 
Channel a acquires a phase relative to channel b. Measure- 
ments are performed on the output state yielding an estimated 
value, ipest, of the phase (p. The beam splitters symbolize pho- 
ton losses. 



measurement procedure that allows one to achieve the 
optimal precision, in terms of a positive operator-valued 
measure (POVM), and the possibility of using states with 
indefinite photon number or distinguishable photons. We 
show that neither of these generalizations improves the 
estimation precision, and consequently the state with def- 
inite photon number and indistinguishable photons are 
optimal. 

The paper is organized as follows. In Sec. [Til a gen- 
eral scheme of quantum phase estimation is presented, 
and the notion of optimality is defined. In Sec. IIIII we 
introduce the quantum Fisher information and discuss 
its most important properties. In Sec. II VI we derive an 
explicit formula for the quantum Fisher information and 
prove that it is a concave function of input state param- 
eters. In Sec. E] we discuss the structure of the opti- 
mal states for interferometry and compare them to al- 
ternative strategies and states. In Sec. I VII we discuss a 
measurement with which it is possible to achieve optimal 
precision. Finally, in Sec. IVIII we discuss possible gener- 
alizations of the considered quantum states, particularly 
states with indefinite photon number and the case when 
photons are distinguishable. 



II. PHASE ESTIMATION 

We consider a general interferometer with two arms 
as shown in Fig. [TJ A pure input state {ipin) is fed into 
the interferometer and acquires a phase ip in the channel 
a relative to the channel b. Both channels, or "arms" , 
of the interferometer are subject to photon losses which 
can be modelled by fictitious beam splitters inserted at 
arbitrary locations in both channels. The output of the 
interferometer therefore needs to be described in general 
by a mixed state p{ip). A measurement, represented by 
a positive operator valued measure (POVM) {11^} which 
defines a probability distribution for the measurement 
outcomes. 



p{i\(p) = TT[n^p{ip)] 



(1) 



ing an estimator that assigns to a particular measurement 
result i an estimated value (^est(*)- We aim to estimate 
the phase (p as precisely as possible. 

Two more elements need to specified in order to make 
the problem of finding the optimal phase estimation 
strategy well defined: an a priori knowledge on the phase 
distribution p{ip) and a cost function C{ip,ipest) which 
can be seen as a measure for the uncertainty in the es- 
timated phase. For phase estimation, the optimal choice 
of the input state, measurement and estimator is the one 
that minimizes the average cost function 



(C) = / d(pp{(p)'^p{i\(p)C{(p,(pest{i))- 



(2) 



At this point two different approaches are most often 
pursued. In the global approach one assumes initial ig- 
norance about the actual value of ip, which corresponds 
to the choice p{(p) — 1/2-k. The solution of the prob- 
lem then yields an estimation strategy which performs 
equally well irrespectively of the actual value of the esti- 
mated phase. In the local approach, on the other hand, 
the assumption is that the value of the actual phase lies 
in the vicinity of a known phase lpq. More precisely, the a 
priori probability is chosen to be p{ip) = S{ip— (pn), while 
the estimator is required to be locally unbiased [3l| 



dip V 



= 1. 



(3) 



is subsequently performed on the output state p{ip). An 
estimated value of the true phase ip is obtained by apply- 



The above condition is equivalent to a statement that the 
estimator will on average yield the true value of ip up to 
the first order in {ip — ipa). Notice that without local un- 
biasedness the estimation problem would be trivial (and 
also useless) since in order to minimize (C) in Eq. 
with p{ip) = S{ip — ipo), the optimal choice for the esti- 
mator would be simply ipestii) = <^0j and the choice of 
the measurement would be irrelevant. 

The local approach is useful when we are interested in 
small deviations of the phase from a known one. A sig- 
nificant advantage of the local approach over the global 
one is that for a natural choice of a quadratic cost func- 
tion, there exist explicit lower bounds on (C) based on 
the Fisher information (see Sec. IIIip . In many practical 
situations these bounds are tight and the optimization 
over the measurement and the estimator can be avoided. 
Moreover, the local approach may also be useful in situ- 
ations when there is no a priori knowledge on the phase 
ip. If many copies of a state are given, one can first per- 
form a rough measurement (even not optimal) on a small 
fraction of copies in order to narrow down the range of po- 
tential values of ip so that they lie in a vicinity of a known 
phase ipo, and then perform the optimal estimation us- 
ing the local approach. This strategy will yield a high 
accuracy estimation for the global approach, without the 
need of optimizing the measurement and the estimator, 
since the rough measurement performed on small frac- 
tion of copies, even if not optimal, will not significantly 
influence the final accuracy j22l |. 
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III. FISHER INFORMATION 

In what follows we take the local approach. Since 
in this case we deal with small deviations of esti- 
mated phases from the true one, it is natural to choose 
C((/3, — {^p — ^p'Y as the cost function. Our goal is to 
find the optimal state, measurement and locally unbiased 
estimator minimizing the expression 



(4) 



To minimize the standard deviation 8ip given by the 
above formula we use a n up per bound on based on the 
Fisher information [lllli]. For a given POVM {HJ and 
state p{}p) defining the probabilities the Cramer- 

Rao inequality bounds the variance that can be obtained 
using any locally unbiased estimator. 



where the Fisher information F is given by 

2 



1 



F = y 



dp{i\ip) 



(5) 



(6) 



If an experiment is repeated v times the bound reads 

1 



vF 



(7) 



For large v the Cramer-Rao bound is asymptotically 
achieved by the maximum likelihood estimator [18, 20, 

Optimization over the measurements yields the quan- 
tum Cramer-Rao bound fisl, HI, ^ [U 



> > 



(8) 



F - Fq' 

where the quantum Fisher information Fq is given by 

Fq = Tr[p{^)A\ (9) 

The Hermitian operator A is called the "symmetric log- 
arithmic derivative" (SLD) and is implicitly defined via 
the relation 



In the eigenbasis of A is given by 

2 



Pi + Pj 



(10) 



(11) 



where p'{ip) — and the pi are the eigenvalues of 

p{ip) (whenever pi + pj = we set {A)ij =0). It has 
been shown that a measurement saturating the quantum 



Cramer-Rao bound exists and is given by a projective 
measurement on the eigenbasis of A [13, [2l| . 

For the sake of completeness we state some important 
properties of Fg (see e.g. [is', '24^): 

(i) Let p{ip), (^{ip) be two density matrices supported 
on orthogonal subspaces, Sp(^,^)J-Sa-(ip), which do not 
cease to be orthogonal for an infinitesimal change of ip, 
i.e. (5p(<^) lJ5p/(^))_L(5c,(^) lJ5c,/(<^)), then Fq is linear 
on the direct sum 



FQ[ppi^)®{l 



-p)a{ip)] 

= pFQ[p{f)] 



il~p)FQ[aiip)]. 



(12) 



(ii) Fq is convex 



Fq [pp{^p) 



(l-p)a(^)] 

<pFq[p{^)] 



{l-p)FQ[a{ip)] 



(13) 



(iii) For pure states p{ip) = \ip{ip)){'ilj{ip)\, Fq reads 



Fn 



4[(V/(^)|^'(^))-|(^'(^)|V'(¥')>I' 



(14) 



where \^'[ip)) — d\ip{ip)) / dip. 

Property (i) is due to the fact that the SLD for 
ppi^p) ® {l—p)a{<p) is a direct product of SLDs for pp{p) 
and (1 — p)a{p), respectively. Property (ii) is a conse- 
quence of the fact that by (i) the right hand side of (ii) 
can be viewed as the quantum Fisher information of the 
state pp{p)(g)\0){0\ + il-p)a{ip) ® |1)(1|, where |0), |1) 
are orthogonal ancillary states, while the left hand side 
is Fq of the state after tracing out the ancillary sys- 
tem. Furthermore, Fq is non-increasing under stochastic 
operations !25| (tracing out the ancilla is an example). 
Property (iii) is a consequence of Eqs. (|8llip . since for a 
pure state 

A = 2(|V'(^))(^'(^)| + \^P'{ip))mip)\). (15) 

The measurement saturating the quantum Cramer-Rao 
bound in this case is a von Neumann measurement pro- 
jecting on any orthonormal basis containing two vectors 



1 



|e±) = -^(|V'(¥'))±l^'^(^))) 



(16) 



where 



i^'^M) = j^iWi^)) - w^wipmip))) (17) 

is the normalized vector orthogonal to \tp{(p)) lying in the 
space spanned by and \'ip'{(p)). 



IV. INTERFEROMETRY WITH LOSSES 

Assuming that we have N photons at our disposal, we 
aim to find the input state that allows performing phase 
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estimation with the best precision possible, i.e. yielding 
the highest value of the quantum Fisher information Fq . 
In particular we consider the most general pure two-mode 
input state with definite photon number N, 

N 

\^p^n)=Y.'^k\k,N-k), (18) 

fe=0 

where \k,N — k) abbreviates the Fock state \k)a\N — k)b- 
This class of states includes the NOON state which, in 
the absence of losses, leads to Heisenberg limited preci- 
sion, but is very fragile in the presence of noise. We are 
therefore looking for states which lead possibly to a lower 
precision than the NOON state, but which are more ro- 
bust with respect to photon losses. Moreover, although 
states of the form ((T8)) seem to be a restriction, we show 
in Sec. I VIII that our treatment effectively includes states 
with indefinite photon number. 

In the following subsections we show how states of the 
form (|18p are influenced by photon losses, calculate its 
quantum Fisher information and show that the latter can 
be maximized (thus minimizing Sip) by means of convex 
optimization methods. 

A. Noise model 

Losses are modeled by fictitious beam splitters of trans- 
missivity rja, r]b in channels a and b respectively, and 
cause a Fock state |fc, iV — fc) to evolve into 

k N-k 

\k,N^k) ^Y.Y1 y/Bl^Jk-L,N-k-lb)<»\la,lb), 

l„=0 /5=0 

(19) 

where \la, h) represents the state of two ancillary modes 
carrying la and It photons lost from modes a and b re- 
spectively, while 

(20) 

Including the phase accumulation \k,N — k) i— > 
e^'^'^lkjN — fc) and tracing out the ancillary modes re- 
sults in the output density matrix 

N N-la 

pi^) = E E pu.\^iM)i^iMi (21) 

1^=0 h=0 

where 

|eu.(^)) = -= E ake^'^'jBtiJk-L^N-k-h) 

(22) 

is the conditional pure state corresponding to the event 
when la and lb photons are lost in modes a and b respec- 
tively, and pi^i^ is the normalization factor corresponding 
to the probability of that event. 



Equivalently, the loss process can be described by a 
master equation for two independently damped harmonic 
oscillators with loss rates ja,b — |ln?ya,b|/^, where t is 
time, the solution of which is given by 

oo 

p = Kl^aK^.,bP^^Kl bK\a (23) 

fc,i=0 

with Kraus operators 

i^i,a= (l-??a)m''V/%/l!, (24) 

where a is the annihilation operator for mode a, and anal- 
ogously for mode b. This state acquires a phase through 
the transformation p(Lp) — q-^v^'' a. p^-fpa-'' a _ js^Q^jee that 
thanks to the relation 

^ig-«Vata ^ ^-iva^a^l^Opl (25) 

we can commute the phase operator with the Kraus op- 
erators since the phase terms e~'^f^ cancels out. It is 
therefore irrelevant if photons are lost before, during or 
after channel a acquires its relative phase with respect to 
b. 



B. Calculating the Fisher information 

Using Eq. (|2T|) for the output state, one can calculate 
Fq with the help of Eqs. ^ and pT]) . This requires 
diagonalization of p{ip) which can be carried out in the 
case of one-arm losses. In the more general case of losses 
in both arms an analytic calculation of Fq turns out to 
be infeasible. Nevertheless, we are able to determine an 
upper bound to the quantum Fisher information which, 
although not strictly tight for general input states, is very 
close to Fq for the states we consider in Sec. |Vl 

1. Losses in one arm 

We consider first the case rja = fj, ijb — ^, i-e. when 
losses are present in only one arm. As can be seen from 
Eq. (PT|) . in this case only states l^ij^if)) with h ~ 
contribute to p{(p). Moreover, we have {^10(^)1^1' oi'p)) = 
Sw, hence we can write the output state as a direct sum 

N 

P(^) -0Pio|6o(^))(6o(¥')l- (26) 

1=0 

Making use of Eqs. and we get a formula for Fq 
with explicit dependence on the input state parameters 

Q!fc, 

fq^a Efc^-^"E I (27) 

yk=0 1=0 l^k=l -^k^io J 
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where Xfe = |afcp. In a more compact way the above 
formula can be rewritten as 



(28) 



where a; is a vector containing variables Xk, while the 
elements of the vector b^'"'*'-' and the matrix RC"''') are 
given by 



if la<k<N-lb 

otherwise 



(30) 



2. Losses in two arms 



If losses are present in both arms, then, in the most 
general case, all \S,i^i^{'f)) contribute to pif). States with 
different total number of lost photons, I = la + hi are still 
orthogonal. Using Eq. (fT2)) we can therefore write 



JV 



1=0 



X! Plal-la \S.lal-la if)) (6ai-/a if) I 



.la=0 



, (31) 



where Fq[-] denotes the quantum Fisher information of 
the state in brackets. Notice that states with the same I 
are not necessarily orthogonal. Consequently, the calcu- 
lation of the above expression requires solving an eigen- 
value problem which is not feasible analytically. Never- 
theless, using the convexity of Fq, Eq. (|13p . we obtain a 
bound 

Fq<Fq^J2 E PuM\iiM){iiM\], (32) 

1^=0 1^=0 

which can be calculated explicitly using Eq. (fH|) . It is 
given by 



N 



N N-l 



k=0 



m=0 Z^fe=i -^fe-^ir 



(33) 



or in a compact form. 



(34) 



;„=0 ii,=o 



When losses are present only in one arm (see previous 
paragraph) Fq = Fq. When losses are present in both 
arms, however, the bound is not always tight. The dif- 
ference Fq — Fq originates from the non-orthogonality 
of IC/a't ('/')) foi' a fixed I = la + lb, and physically cor- 
responds to lack of knowledge about how many photons 
were lost from a particular mode. If this knowledge is 
not relevant then Fq = Fq. This happens, e.g. in the 
case of the NOON state ao|^,0) -I- a^vlO, TV), when loss 
of even a single photon renders the output states useless 
for phase estimation, hence the knowledge of which mode 
the photons were lost does not influence the value of Fq . 



C. Concavity of Fisher information 



Even with the explicit formulae for Fq and Fq given 
by Eqs. (p8| and ((34| it is in general not possible to find 
an analytic solution for the optimal state, i.e. values 
Xk > 0, X^fc^fc — 1 that maximize Fq (or Fq). How- 
ever, we prove below that Fq is a concave function of the 
{xk}- Consequently, the problem amounts to the maxi- 
mization of a concave function Fq on a convex set. This 
allows for a feasible numerical constrained optimization 
using interior-point method routines (e.g. implemented 
in Mathematica 6.0), and more importantly, any local 
maximum found is automatically the global maximum. 

We prove concavity by showing that the Hessian Hij = 

is negative semidefinite, i.e. for every vector y 

we have y^Hy < 0. Using Eq. (f34|) the Hessian Hij 
{ij G {0, . . . ,iV}) reads 



_ ^ ^» bf-'^h^^-'^Hx^R^'^'^^x) - a;^b('«'^)[6f°'''^(R('"'^)a;)j + bf-''\R^'^'^^x),] + x^b^'^'^^R^'-'^^.x^b^'^'^ 

ki=0 lt=0 ^ ' 

(35) 

Since the denominator is always positive, it is sufficient to prove that for every 1^, lb and for every vector y we have 
(yT^(i„i.))2(^Tj^(ij,)^) - (a;^6('»''))(y^6('"''))(y^R('-'^)a; + a;^R('-'^)y) + (a;^b('"'^')2y^R('»'^)y < 0. (36) 

I 

Introducing K,j ^ (i - j)^ , Wi = y^fef (aj^b*'"'"^) - x^fof "'"^(y'^b^'"'"^) the above condition can be written 
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equivalently as 



w^Kw < 0. 



(37) 



Since Wi = 0, it is sufficient to prove that K is neg- 
ative semi-definite on the set of vectors with coefficients 
summing up to 0. To this end we define vectors e^*) 



{i = f , . . . , N), where e 



(0 



f, 



-1, e 



(0 



(for 



j 7^ 0, j 7^ i), which span the space of aU vectors with co- 
efficients summing up to zero. Writing w = Afi^*'* 
and noticing that e*^*^"^Ke*^-'^ = —2ij we arrive at 

w'^Kw = -2 ^ (3^(3jij = ~2 I ^ < 0, (38) 

i,j \ i / 

which proves that Fq is a concave function of the {xk}- 



V. OPTIMAL STATES FOR PHASE 
ESTIMATION 

In this section we discuss optimization results based 
on Eqs. (^7)) and derived in the previous sections. 
The quantity we analyze is S(pmin = 1/ \/Fq correspond- 
ing to the best possible precision for a fixed number of 
measurements i/ (cf. Eq. ([8])). The only exception to this 
definition is in the case of the optimal state for losses in 

both arms where we set Sipmin = (see Sec. IV Bt . 

We compare the optimal precision to the precision which 
is obtainable using various alternative states and strate- 
gies. In particular, we define the standard interferometric 
limit (SIL) corresponding to the precision which can be 
achieved in a classical reference experiment. This serves 
as a benchmark by which we can judge the advantage of 
using quantum states of light over classical ones. For a 
given photon number N the SIL is the precision of phase 
estimation using a Mach-Zehnder interferometer in which 
one arm gathers a phase f, fed at one input port with a 
coherent state |a), where |ap = N, and the vacuum at 
the other port. The reflectivity of the two beam splitters 
can be adjusted to achieve the best precision, while the 
measurement consists of photon counting at the two out- 
put ports. Without any additional reference beams the 
input coherent state should be regarded as a state with 
unknown phase, and effectively described as a mixture 
of Fock states with Poissonian statistics (see [1^, and 
the discussion in Sec. IVIlXj) . By Eq. the quantum 
Fisher information will be a weighted sum of quantum 
Fisher information calculated for each input Fock state 
after it passes through the first beam splitter. Taking the 
optimal value of transmissivity of the first beam splitter 
such that the ratio of the intensities in the arms a and b is 
\/ Vb /Va the final uncertainty in phase estimation achiev- 
able with this classical strategy is given by 



dtp. 



mm, SIL 



(39) 



Note that the SIL scales in the same way as the SQL, 
particularly for equal losses in both arms we obtain 
SipsiL = where t] = rja = rjb. 

As a reference let us also calculate the minimum un- 
certainty achievable using N photons prepared in an un- 
balanced NOON state: 



/i^^lOTV). 

The quantum Cramer-Rao bound yields: 



'mm, NOON 



N/2 , N/2 
N/2 N/2' 



where the optimal amplitudes are given by 



N/2 



N/2 



N/2 



(40) 



(41) 



(42) 



and xn = 1 — xq. Note that putting A'' = 1 in Eq. ([4T|) 
coincides with Eq. ([59)1 . This implies that a coherent 
state performs equally well in phase estimation as the 
equivalent number of single photons sent one- by-one. In 
the following we will concentrate on the two important 
scenarios of losses only in channel a and equal losses in 
both channels. 



A. Losses in one arm 

The scenario of losses in only one arm of the interfer- 
ometer, i.e. rja — rj and rji, = 1, is relevant, for example, if 
losses are induced by a sample itself. Figure [2] shows the 
parameters xu of the optimal state for iV = 10 photons 
as a function of r\. The corresponding precision is shown 
in Fig. m The shaded, grey area in this figure is bounded 
by the SIL and the Heisenberg limit 1/iV, i.e. whenever 
a line is in this region there is an improvement over clas- 
sical interferometry. The precision obtainable with the 
NOON state is worse than the SIL except for relatively 
low losses. 

For transmissivities exceeding a certain threshold r\ > 
f\ the optimal state is equal to the NOON state (140]) . Note 
that for ry = 1 the NOON state is "balanced", i.e. = 
x-N = 1/2. We will give an estimate for the threshold 
fj at the end of this subsection. For t] < fj the optimal 
state consists of more components. As seen in Fig.[3]over 
a wide range we can achieve almost the same precision 
by using two-component states of the form 



/^|fc,iV- 



(43) 



where we numerically optimize k and the amplitude Xk- 
With decreasing r] the optimal choice for k increases. 
This state is more robust in the presence of losses than 
a NOON state, since a loss of up to k photons in the first 
mode does not destroy the superposition, while for the 
NOON state a loss of even a single photon renders the 
state useless for phase estimation. On the other hand. 
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FIG. 2: (color online) . Parameters of the optimal state: j^/") = 
X]i=o N — i), for phase estimation with A'^ = 10 photons 

for the case of losses in one only arm, i.e. rja — rj, rjt = 1. 
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FIG. 3: (color online). The minimal phase uncertainty 
achieved with various A'^ = 10 photon states for losses in one 
arm, i.e. i]a — ri, rjb = 1. Red, dashed line: NOON state; blue, 
solid line: optimal state (cf. Fig. [5}; bright green, dotted line: 
optimal two-component state; dark green, dashed-dotted line: 
NOON chopping strategy. The shaded area is bounded by the 
SIL achievable with classical states and the Heisenberg limit 
1/iV. 



a large N — k increases the sensitivity of the state with 
respect to an induced phase in channel a. The properties 
of large k and large N — k are therefore competing and 
the optimal result represents a trade-off between phase 
sensitivity and robustness. 

An alternative strategy that leads to an improvement 
over the SIL, but uses states of a simpler structure than 
the optimal state is the NOON "chopping" strategy which 
was introduced in Ref. [l^. Instead of a single NOON 
state, we use the same number of photons, but send them 
successively in N/ n smaller portions using n-photon nOOn 
states. Repeating an experiment N/n times corresponds 
to an N/n fold increase in Fisher information, as shown 
in Eq. ([7]). Using Eq. ([ill) we find that the Fisher infor- 



mation for the NOON chopping strategy reads 



i'O — — 



n (l + 77"/2)2- 



(44) 



Treating both n and N/n as real numbers with n G [1 , -/V] , 
maximization of this expression over n yields 



N 
N 



In 71 

In r)a 



;?7 < ?7o ~ 0.228 
■,Vo <V<Vo 



(45) 



where the optimal choice of n for the three regimes indi- 
cated above is n = 1, n = 1.478/| \nr]\ (i.e. the solution 



of 1 



n/2 



-|- nln?7 = 0) and n = N. We see that for 
very small transmissivities t] < rjQ the precision is equal 
to the SIL, specified in Eq. (|39p. For higher transmissiv- 
ities the chopping strategy beats the SIL, although only 
by a constant factor rather than in terms of scaling. For 
very high transmissivities rj > t]q^^ the best strategy is 
to use un-chopped NOON states. We note that a similar 
strategy for multipartite qubit states has been devised 
in [12]. 

As mentioned above the NOON state ceases to be op- 
timal below a threshold fj in which case the coefficient 
xi obtains a nonzero value. By adding an infinitesimal 
change 6xi to xi at the expense of xo,xn (which we as- 
sume are kept in the proportion which is optimal in the 
absence of .ti), we can calculate the corresponding change 
in the quantum Fisher information. The change of xi by 
6xi increases Fq by ^^Sxi, but at the same time due 
to decreasing weight of the other coefRcients decreases it 
by FqSxi. On the whole the change of Fq reads: 



SFc 



dFQ 
dxi 



Fq Sxi. (46) 



Substituting Fq from Eq. ^ Sit xa ^ ri^l'^/{\ ^ 
xm = 1 — xq, Xi — {i ^ 2, . . . , N — 1) and calculating 
5Fq we determine the value 77 below which SFq/Sxi is 
positive. This implies that for ij < fj an increase in xi 
will increase Fq. Writing 6Fq/6xi explicitly we get: 



SFq 



[(1 - 27V)ry^ - 2N{N-l)r] 



N/2+1 



Sxi (1 + 77^/2)2 ' 

+ 2{N - 1)^77^/2 _ _ 2)?7 +{N- 1)2]. (47) 

The roots of the expression in square brackets can be 
found numerically. There is one real root in the interval 
[0, 1] which corresponds to the threshold fj. The roots 
behave in very good approximation like fj — a^^/^ where 
a « 2.61 which is obtained by a fit to the roots between 
N = b and N = 100. For example, for A'' = 10 we find 
that the threshold at which the NOON state ceases to be 
optimal is at 77 = 0.91, which agrees with the numerical 
results obtained before (see Fig. [2). 
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Parameters of the optimal state: — 
i,N — i), for phase estimation with A*' = 10 photons 



for the case of equal losses in both arms, i.e 
Due to this symmetry we have Xk = XN-k- 



m = V- 



B. Equal losses in both arms 

In order to determine the optimal state in the case 
of equal losses in both arms of the interferometer, i.e. 
r]a — rjb = 1], we use Fq as given by Eq. as a ba- 
sis for our numerical optimization. It simplifies maxi- 
mization drastically since it is a concave function of the 
{xk}. As pointed out in Sec. lIVB"2l Fq, is only an upper 
bound to the "true" quantum Fisher information Fq and 
therefore cannot be reached by any measurement strat- 
egy if Fq is strictly smaller than Fq. However, using Fq 
can be still very useful: We maximize Fq and use the 
resulting input state (say, \ipin)) to calculate the corre- 
sponding FQ^ipin)). On the one hand, the true maxi- 
mum of FQ^ipin)) (corresponding to the truly optimal 
state, I "04™)) is certainly greater or equal to FQQipin)). 
On the other hand we have FgdV'in)) < -PgdV'm)) due 
to convexity [see Eq. ([5^ ]. But since -FgdV'm)) is the 
global maximum we have fgd^'in)) < -pQdV'm))- Conse- 
quently, the true maximum must lie in between -FgdV'm)) 
and FgdV^in)). Hence, if iFgdV^^n)) - i^Q(|-0j„))| is small 
on a scale which is for this problem naturally given by the 
difference between the SIL and the Heisenberg limit, we 
obtain a very good approximation of the true maximal 
value of the quantum Fisher information. We showed nu- 
merically that for the examples discussed in this section 
this difference never exceeds 0.4%. 

The result of optimization is shown in Fig. [5] for N = 
10. Due to the symmetry of the problem, rja — rji,, the 
solution has the property Xk = x^^k- For small losses 
the optimal state is the balanced NOON state, i.e. xq — 
xn = 1/2. However, below a threshold fj all other terms 
Xk obtain non-vanishing values. With decreasing 77, xq 
and Xn get smaller, and eventually the coefficients Xk 
with k « N/2 are dominant. Interestingly, this structure 
is qualitatively similar to a two-mode generalization of 
the state which is optimal for one-mode lossless phase 



FIG. 5: The minimal phase uncertainty achieved with var- 
ious A'^ = 10 photon states for equal losses in both arms, 
i.e. ria = rit = rj. Red, dashed line: NOON state; blue, solid 
line: optimal state (cf. Fig. |4)|; bright green, dotted line: 
the state optimal for lossless phase estimation in the global 
approach [see Eq. (|48p ]: dark green, dashed-dotted line: the 
NOON chopping strategy. The shaded area is bounded by the 
SIL achievable with classical states and the Heisenberg limit 
1/iV. 



estimation in the global approach [27| 

N 

IV^™) = ^2/{N + 2) ^sin[7r(fe + l)/(iV + 2)]\k, N-k). 

k=0 

It was noted in Ref. ^28] that these states perform quite 
well in the presence of loss, and hence could be useful for 
the problem we consider. 

The threshold fj can be obtained in a similar fashion 
as in the previous subsection. By adding an infinitesimal 
change Sxk to a component Xk other than xq and xn we 
can derive a polynomial which has a root in [0, 1] defining 
fj{N). Since the polynomial is rather cumbersome we do 
not quote it here, however, the relevant root is again 
given in very good approximation by 77 = a~^/^ where 
a « 2.24 is obtained by a fit to the roots between N = 2 
and = 100. For example for = 10 (corresponding 
to Figs. Hand O we get fj = 0.92. 

Analogous to the case of losses in oidy one arm we can 
examine a NOON chopping strategy by sending "smaller" 
nOOn states through the interferometer. The correspond- 
ing precision reads 



c| 111 ?7| 

N 



; ?7 < e ^ 



; e ^ < 7/ < e 



(49) 



;?7 > e" 



FigurelH shows the estimation precision for the optimal 

state (blue line; here we use Sipmin = 1/ \J Fq), the NOON 
state (red line), the NOON chopping strategy (dark green 
line), and the state (|48p (bright green line). In the latter 
case the precision is indeed quite high in the regime of 
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high losses, outperforming the NOON chopping strategy 
but for smaU fosses it is significantly worse than both the 
optimal and the NOON chopping strategies [3^ . 

VI. THE OPTIMAL MEASUREMENT 

As discussed in Sec. IIIH a measurement saturating the 
quantum Cramer-Rao bound always exists. For the cases 
when the states {(.lahi'P)) introduced in Sec. IIVI are or- 
thogonal for different I ~ la + h (i.e. when Fq — Fq) 
we can derive an explicit POVM saturating the bound. 
To this end we first assume that we perform a measure- 
ment determining the total number of photons lost, /, 
which projects the system onto a particular pure state 
\^iaib{'-P)) |S3l. Subsequently, we perform the measure- 
ment saturating the Cramer-Rao bound on pure states 
(see Sec. IIII[) . i.e. a projection on a basis containing the 
vectors 

where (na)i„ib, (Att-q)^ are the mean photon number 
and the variance of the photon number in mode a for the 
state l^.i^ibi'p)) ■ The above POVM depends, in general, 
on the phase (p. 

There are cases in which a single POVM saturat- 
ing the Cramer-Rao bound for all values of (p can be 
found. In a single photon case {N — 1), states of the 
form l/^/2(|0,l) + e^^\l,0)), allow for a 93-independent 
POVM, whereas states with unbalanced weights of |01), 
and 1 10) terms lead to optimal POVMs which are nec- 
essarily (^-dependent. Moreover, for an N photon state 
IV') = J2k=o^k\k,N — k), a single, (^-independent POVM 
saturating the Cramer-Rao bound can be found pro- 
vided that the state enjoys the path-symmetry property: 
ttfc = c(%-k^^^^ where x is a fixed phase the same for 
all components [1^. In general, due to losses the con- 
ditional pure states that appear in our problem lack the 
above mentioned symmetry and consequently the opti- 
mal POVM varies with the change of (p. 

VII. MORE GENERAL INPUT STATES 

One could argue that the N photon state given by 
Eq. docs not represent the most general way of us- 
ing N photons to determine the relative phase in an in- 
terferometer: Instead of considering a state with a defi- 
nite photon number N, one could consider a state which 
is a superposition of terms with different total photon 
number, and only imply a constraint that the mean pho- 
ton number is N. Furthermore, by using states of the 



form (jlSp it is assumed that all photons in an arm of the 
interferometer are indistinguishable, i.e. occupy a com- 
mon spatio-temporal mode. This greatly limits the class 
of considered states since the Hilbert space dimension of 
states sent into the interferometer would be of the order 
2^ if the photons were distinguishable, while for indistin- 
guishable photons it is iV -I- 1. In the next two subsections 
we show that none of these generalizations improves the 
phase estimation precision, and consequently it is suffi- 
cient to consider the states of the form given by Eq. (jlSp . 

A. States with indefinite photon number 

We emphasize that in this paper we consider closed 
systems, i.e. there are no additional reference beams, nei- 
ther classical nor quantum, since any additional beam 
would contain photons and should therefore be explicitly 
taken into account for the determination of the required 
resources. 

Consider a superposition of input states of the inter- 
ferometer with different, but definite photon numbers n, 

\^.n)=f:Pn\4% (51) 

(50) 

where the mean total photon number is = 
X]^o 1/5" P"-- Under free evolution a term with n pho- 
tons acquires a phase e~™'^*, i.e. terms with different 
n evolve with different frequencies. In the absence of an 
additional reference beam which allows for clock synchro- 
nization between the sender and the receiver the relative 
phases between terms with different n become unobserv- 
able and the state given by Eq. (|?T|l is physically equiv- 
alent to a mixture |26l . [30| , 

00 

P^n = Y.\f^M:J){4:^\■ (52) 

Moreover, since the quantum Fisher information is con- 
vex (see Sec. IIIip. we have 

00 

fq(p„) < m'FQ{\4:^){4:^\). (53) 

Consequently it is always better to send a state with a 

fixed number of photons IV"!-"^)? with probability |/3nP, 
rather than to use a superposition (which is effectively a 
mixture). The analysis can thus be restricted to states 
with definite photon number without compromising op- 
timality. 

B. Distinguishability of photons 

In this subsection we consider general N photon input 
states where the photons are distinguishable which is the 
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case, e.g., if they are sent in different time bins. A state 
of this type can be written as 

1 

|V'™) = 5I"fel^)' (54) 

k=0 

where k represents a binary sequence of length N, and 
= {0, . . . , 0}, 1 = {1, . . . , 1}. Summation from a se- 
quence to a sequence should be understood using partial 
ordering between sequences, i.e. k < k' ki < k'^ Vi. 
Therefore the summation is performed over all binary se- 
quences. Furthermore \k) = \ki) (g) • • • ® |fcAr) denotes a 
state of N photons, where ki — 0(1) means that the pho- 
ton sent in the i-th time bin propagates in arm a (b) of 
the interferometer. Mathematically, we deal here with an 
arbitrary 7V-qubit state, which lives in a 2^ dimensional 
space. Notice, that the indistinguishable case is recov- 
ered once we consider only the fully symmetric (bosonic) 
subspace of N qubits. 

Fq can be calculated in a similar way as in Sec. IIVBI 
and reads 

Fq=A Efe"^-E E y^i-^., h 

(55) 

where boldface symbols denote binary strings of length 
A'', Xk = |afcP, k denotes the number of Is in the se- 
quence fe, and 

Bti, - (1 - VaY-vt'-il - (56) 

Again Fq — Fq provided that knowledge from which 
mode a photon was lost is not relevant or self- 
evident, e.g. when rji, = 1. Let us now sym- 
metrize the state of N photons. This corresponds 
to replacing Xk with 1/A^! X^o- ^o-(fe)' where the sum- 
mation is performed over all permutations of an A^ 
element set. Notice that the first term within the 
parenthesis in Eq. (|55p is not affected by symmetriza- 
tion. Let us denote aj^j^ = ^k=\'' ^^kB^ and 
^lah — Efe^j'' ^kBf I . Then the subtracted term in 

2 

Eq. dSS]) reads EL=o Ei,=o ^rflf"- Performing sym- 

metrization corresponds to replacing ai^i^ and with 

l/^!EaMia),^(ib) and l/N\Y,^ K{i^),a{h) respectively. 
Since bi^ii^ are positive, the following inequality holds 



which shows that symmetrizing can only decrease the 
subtracted fraction, hence increase the Fisher informa- 
tion. This proves that the optimal states are symmetric 
states living in the bosonic subspace. Consequently, us- 
ing indistinguishable photons is sufficient to obtain the 
optimal Fisher information. 

A physical intuition behind the above derivation is the 
following. Notice that if there are no losses, the optimal 
A^ qubit state has the form l/\/2(|0) -I- |1)) - it is the 
NOON state which lives in the fully symmetric subspace. 
In this case there is no advantage in using distinguishable 
photons. If there are losses, things get even worse for 
distinguishable photons. There is still no advantage in 
terms of phase sensitivity, yet when a photon is lost, the 
knowledge of which photon was lost additionally harms 
the quantum superposition. Hence, it is optimal to use 
states from a fully symmetric subspace, where it is not 
possible to tell which photon was actually lost. 



VIII. SUMMARY 

We have analyzed the optimal way of using A^-photon 
states for phase interferomctry in the presence of losses. 
We have derived an explicit formula for the quantum 
Fisher information in the case when losses are present 
in one arm of the interferometer, and have provided a 
useful bound for the case of losses in both arms. Using 
the quantum Fisher information as a figure of merit, we 
have found the optimal states, investigated their advan- 
tage over various quantum and classical strategies. The 
optimal measurement saturating the Cramer-Rao bound 
has been presented, and it has been proven that in gen- 
eral it varies depending on the phase in the vicinity of 
which we perform estimation. 

A close inspection of the properties of the quantum 
Fisher information showed that it is optimal to use a 
quantum state with a definite number of indistinguish- 
able photons. Neither allowing superpositions of different 
total photon-number terms, nor making photons distin- 
guishable can improve estimation precision. 
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